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ABSTRACT 


A  three  dimensional  finite  element  code  is  written  for 
truss  analysis  and  design.  Trusses  may  be  designed  for  mini¬ 
mum  weight  subject  to  constraints  on:  member  stresses,  Euler 
buckling,  joint  displacements  and  system  natural  frequency. 

The  optimum  configuration  may  be  found  in  addition  to  optimi¬ 
zation  with  respect  to  member  sizes. 

The  finite  element  code  may  be  used  as  a  stand  alone 
analysis  tool  or  may  be  coupled  to  an  optimizer  of  the  user's 
choice.  The  finite  element  displacement  method  of  analysis 
is  used  for  static  analysis  and  eigenvalues  are  calculated 
using  the  subspace  iteration  technique. 

A  major  portion  of  this  work  is  the  software  (user  guide) 
which  is  presented  in  detail  with  examples  and  results.  Ex¬ 
planation  of  how  this  code  is  coupled  to  an  optimizer  is 
given. 

Design  variables  are  member  area  sizes,  joint  coordinates, 
or  both.  Examples  are  presented  to  demonstrate  the  method. 
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I.  INTRODUCTION 


Design  optimization  requires  the  minimization  or  maximi¬ 
zation  of  some  parameter.  Optimization  of  structures  has  had 
continuing  changes  since  its  development  in  the  early  1960 ' s 
with  an  active  area  of  research  being  elastic  truss  structv  s. 
The  main  goal  is  to  design  structural  systems  that  efficien  v 
perform  specified  purposes .  For  any  design  to  be  acceptabl 
it  must  satisfy  a  variety  of  physical,  aesthetic  and  economic 
constraints . 

Since  most  physical  problems  can  be  modeled  by  some  mathe¬ 
matical  formulation  a  computer  program  can  be  written  to  per¬ 
form  the  necessary  calculations . 

The  purpose  of  this  research  was  to  develop  a  FINITE 
ELEMENT  code  that  could  be  easily  coupled  to  an  optimizer, 
thus  (1)  allowing  testing  and  comparison  of  various  optimizers 
and  (2)  provide  a  useful  design  program  in  its  own  right. 

The  design  problem  considered  in  this  study  is  the  opti¬ 
mization  of  three-dimensional  indeterminate  trusses ,  for 
multiple  static  load  conditions.  The  objective  is  to  mini¬ 
mize  the  weight  of  the  structure  where  the  design  variables 
are  member  sizes  and  joint  coordinates.  Constraints  include 
stress,  displacement,  Euler  buckling,  and  the  natural  fre¬ 
quency  of  the  system. 
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This  document  describes  the  use  and  capabilities  of  the 
finite  element  computer  code  to  be  coupled  to  an  optimizer. 

The  user's  manual  presented  in  Chapter  V  contains  a  simple 
design  example  in  which  the  program  is  coupled  to  the  CONMIN 
optimization  code  [Ref.  l]. 

Additionally,  guidelines  for  coupling  the  code  to  an 
optimizer  of  the  user's  choice  are  presented. 

Several  examples  demonstrating  the  program  under  a  variety 
of  conditions  are  presented.  Conclusions  and  recommendations 
for  future  work  are  given. 
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OPTIMIZATION 


II  . 

A.  INTRODUCTION 

The  main  goal  of  structural  engineering  optimization  is 
to  design  structural  systems  that  efficiently  perform  speci¬ 
fied  purposes.  Selection  of  a  specific  algorithm  must  include 
the  following  considerations:  1)  the  structure  should  be 
analyzed  as  few  times  as  possible,  2)  the  algorithm  should 
minimize  specific  gradient  information  required,  and  3)  it 
should  provide  reasonable  assurance  that  an  optimum  design 
will  be  reached. 

The  next  few  sections  discuss  the  general  formulation  of 
the  algorithm  with  respect  to  the  above  requirements. 

3 .  FORMULATION 

Minimum  weight  design  of  trusses  is  presented  in  the 
general  form  of  a  mathematical  programming  problem  as  follows: 

Minimize  F(X)  (Eq.  1) 

Subject  to: 

gj  (x)  _<  0  j=l,m  (Eq.  2) 

u 

xi  £  xi  1  xi  i=l,n  (Eq.  3) 

where  F(X)  is  the  objective  function,  in  this  case,  weight  of 
the  structure  be  minimized.  X  is  the  vector  of  independent 
design  variables  and  here  contains  the  member  cross-sectional 
areas  as  well  as  the  coordinates  of  the  joints.  The  inequality 
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constraints,  g^(x)  <_  0,j=l,m,  must  be  satisfied  for  the  design 

to  be  accepted  as  feasible.  These  include  limits  on  stress, 

Euler  buckling,  joint  displacements  and  the  first  fundamental 

.  0 

frequency  of  the  structure .  Side  constraints  x^  and  x^  are 
lower  and  upper  bounds  on  the  design  variable.  These  may  be 
treated  as  general  inequality  constraints. 


C.  PARAMETERS 

The  design  variables,  constraints,  and  objective  function 
considered  in  the  design  process  are  discussed  here. 

1 .  Design  Variables 

The  design  variables  the  user  has  available  are  member 
cross-sectional  area,  reciprocal  of  member  cross-sectional 

area,  joint  coordinates,  or  both  member  areas  (or  their  reci-  ; 

! 

procals)  and  coordinates.  The  cross-sectional  areas  are  A^ , 

k=l,NE,  where  NE  is  the  number  of  elements  in  the  truss  i 

I 

structure.  The  joint  coordinates  are  ,  i=l,2,3  j  =  l,  NJ,  j 

where  i  is  the  coordinate  axis,  j  is  the  joint  number  and  NJ  j 

is  the  total  number  of  joints.  Design  variable  linking  is 
allowed  in  both  member  sizing  and  coordinate  design  variables. 

2  .  Objective  Function 

The  objective  function  considered  here  is  weight. 

NE  AL 

Weight  (W)  =  Z  p.  bi  (Eq.  4) 

i  =  l  1 


where  p  is  the  material  density,  A  is  the  cross-sectional  area 
and  is  the  length  of  the  member.  The  truss  may  be  made  up 
of  members  of  differing  materials. 
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3 .  Constraints 

i = 1 , ME  (Eq.  5) 

a .  . 

Stress:  -  1  <_  0  j=l,MLC 

5 


-  1  <  0 


i=l ,NE  (Eq.  6) 

j  =  1 ,MLC 


where  a  is  the  stress  in  member  i  under  load  condition  j .  NE 
is  the  number  of  elements  and  NLC  is  the  number  of  loading 
conditions . 

a  ~  is  the  lower  bound  on  the  stress  (maximum  compressive 
stress)  and  a  +  is  the  upper  bound  on  the  stress.  The  upper 
and  lower  bounds  on  stress  may  be  different  for  each  member, 
but  are  taken  to  be  the  same  for  every  loading  condition. 

4 .  Euler  Buckling 

The  stress  at  which  Euler  buckling  occurs  is  given  by 


a 


-KiAiEi 


i = 1 , NE  (Eq.  7) 


where  the  subscript  i  corresponds  to  the  member  number,  E^ 
is  Young's  modulus  and  is  a  constant  depending  on  the 
cross-sectional  geometry  of  the  member. 

5 .  Displacement 

Displacement  limits  are  imposed  at  prescribed  joints 
to  create  a  constraint  equation  as  follows: 
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jjjJl  _  1  <  0 

uijk 


(Eq.  8) 


-  i  <  o 
+  — 
u. 
ink 


(Eq.  9) 


where  u.^^  and  u^j£  are  lower  and  upper  bounds  on  the  displace¬ 
ment  at  joint  i  in  the  coordinate  direction  j  under  loading 
condition  k. 

6 .  Frequency 

The  first  fundamental  frequency  of  the  structure  is 
required  to  exceed  the  specified  lower  bound,  so  that 


1  -  X/X  <0 


(Eq.  10) 


Reference  2  is  an  excellent  source  for  the  basic 
structural  design  formulation. 


15 


III.  FINITE  ELEMENT  METHOD 

A.  INTRODUCTION 

Several  features  are  desirable  when  the  finite  element 
methods  of  analysis  are  used  in  design  optimization.  First 
the  number  of  analyses  for  the  structure  should  be  kept  to  a 
minimum.  Second,  the  amount  of  gradient  information  required 
during  the  design  process  should  be  reduced  to  shorten  run 
times  and  computer  storage  requirements.  Third,  the  user 
should  be  able  to  specify  only  that  gradient  information 
desired . 

3 .  ANALYS IS 

Initial  problem  formulation  includes  member  sizing,  mate¬ 
rial  properties  (which  may  be  different  for  each  member) ,  a 
set  or  sets  of  external  loads,  and  specified  support  condition. 

The  analysis  for  the  stresses  and  deflections  must  satisfy 
the  conditions  of  equilibrium  of  forces  at  the  nodes  and  com¬ 
patibility  of  deformation.  In  this  analysis  the  weight  of  the 
individual  members  are  not  part  of  the  specified  load  conditions. 

Additionally  for  this  analysis  the  following  assumptions 
are  made.  Trusses  will  be  treated  as  discrete  elements,  and 
each  element  will  be  treated  as  pin-connected  with  loads  and 
reactions  supported  at  the  joints. 
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The  general  method  of  solution  is  as  follows.  The  Dis¬ 
placement  (Stiffness)  method  reference  2  considers  the  joint 
displacement  components  as  the  unknowns,  and  written  in 
matrix  notation  for  the  general  case  is 

Ku  =  P  (Eq.  1) 

where  K  is  the  global  stiffness  matrix,  P  is  the  vector  or 
vectors  of  applied  loads,  and  u  is  the  vector  of  vectors  of 
displacements . 

Once  displacements  at  every  node  are  known,  the  internal 
forces  and  stresses  are  calculated  by  applying  the  appropriate 
force-deflection  relationships. 

When  the  system's  natural  frequency  constraints  are  con¬ 
sidered  the  design  process  requires  the  solution  of  an  eigen- 
problem.  This  solution  will  determine  the  natural  frequencies 
and  normal  modes  of  the  structure.  For  linear  elastic  struc¬ 
tures,  the  finite  element  approach  leads  to  the  following 
equation  of  motion  with  free  vibration  conditions, 

Mu  +  Ku  =  0  (Eq.  2) 

where  M  is  the  global  mass  matrix,  and  u  is  the  linear  accel¬ 
eration  vector.  This  leads  to  an  eigenvalue  problem  of  the  form 
IK—  X  ld=0 ,  and  is  solved  here  by  the  subspace  iteration  method 
reference  2  . 
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C.  GRADIENTS  OF  CONSTRAINTS 

Gradient  computation  of  necessary  functions  in  a  design 
optimization  process  arise  from  the  need  for  derivative  infor¬ 
mation  for  efficient  mathematical  programming.  For  structural 
optimization  the  gradient  of  displacement  with  respect  to  the 
design  variables,  3u/3x,  is  needed  from  which  VG ,  is  calculated 
for  stress,  displacement,  and  buckling. 

Consider  the  finite  element  method  Ku  =  P .  Taking  the 
derivative  of  both  sides  with  respect  to  the  design  variable 
x£ ,  we  have : 

—  CKu]  =  —  (P)  (Eq.  3) 

3xz  =~  3x£ 

Assuming  the  loads  P  are  not  a  function  of  x£  , 

[—  (K)]  (u)  +  K  —  (u)  =  0  (Eq.  4) 

3x*  *  ~  =  3*a 

Finally,  we  arrive  at  (u)  =  -K-1  [— —  (K)]  u  (Eq.  5) 

3x£  ~  =  3x£  s 

where  K~^  is  the  inverse  of  K.  For  efficiency  K-^"  is  not 
actually  calculated. 
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IV.  PROGRAM  FEATURES 


A.  INTRODUCTION 

Each  computer  code  has  its  own  special  features  with  which 
the  user  should  be  familiar  if  the  program  is  to  be  used 
effectively.  The  next  few  sections  will  discuss  an  overview 
of  the  code  and  a  typical  problem  that  might  be  solved  using 
the  program.  This  problem  along  with  other  numerical  examples 
will  be  presented  in  detail  with  results  i~  Chapter  VI. 

The  FINITE  ELEMENT  code  was  written  to  be  used  as  a  stand 
alone  analysis  program  or  an  analysis  code  that  could  easily 
be  coupled  to  an  optimizer  (of  the  user's  choice)  through 
simple  modifications  to  the  main  driver  program. 

With  user  supplied  area  and  coordinates  coupled  with  input 
control  parameters,  the  analysis  mode  will  calculate  the  weight 
of  the  truss  structure.  Design  variables  may  be  chosen  as 
member  areas,  joint  coordinates,  or  both.  Additionally, 
gradient  information  will  be  calculated  with  respect  to  area 
variables,  coordinate  variables  or  both.  Coupled  to  an  opti¬ 
mizer  the  code  will  optimize  the  weight  of  the  structure  and 
print  the  final  optimization  information  and  gradient  vectors 
of  specified  constraints. 

The  following  example  of  a  25-bar  space  tower  presented 
in  Tables  (I-IV)  shows  some  of  the  options  the  code  contains. 
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Table  I 


2 5 -BAR  TOWER  SAMPLE  OUTPUT 


INPUT  CONTROL  PARAMETERS 


INPUT  PAR AM£T£SS  FOR  STRUCTURAL  ANALYSIS  AND  DESIGN  ROUTINE,  "SAD" 
25-3AR  SR4CF  TRUSS  (TEST  CASE  STP ESS7DI S P./3UCKLI NS/FRfcU. ) 

control  parameters 


NUMBER  OF  ELEMENTS,  NS  = 

TOTAL  NUMBER  CF  JOINTS,  NJ  3 

JOINT  CONSTRAINT  VARIABLE,  NCJ  = 
NO.  OF  M  A  T  E  R.  I  AL  TYPES.  NMT 
NO.  OF  LOAD  CCNO . •  NLC 
NC.  CF  EIGENVALUES,  NEIG 
NO.  OF  c IGENVALUES  CALC.,  NE IG1* 
NO.  OF  FIXED  MASSES.  NFMASS 
BUCKLING  CONSTRAINT  10,  NEUBC  ’ 
NO.  OF  3ISPL.  CONST., NOSPLC 
NO.  OF  FRFO.  CONST.,  NFREQ 
LUMPED  MASS  OPTION,  LMASS 
V A°B.  CALC.  CONTROL  ,  IOVCLC 


25 

10 

10 

1 

2 

2 

2 

2 

1 

4 

1 

0 

2 


ACCELERATION  OUE  TO  GRAVITY,  GRAV  _  -  0.38640E*03 

EIGENVALUE  CONVERGENCE  TOLORANCE.  EPSEIG  «  O.IOOOOE-Q3 


JOINT  COORDINATES 
JOINT  X 

1  -0.375006*02 

2  0.37500C*02 

3  -3.3  75006*02 

4  0.375G0E*02 

5  0.375036*02 

A  -0.375006*02 

7  -0. 100006*03 

8  0. 100006*03 

9  0.  100005*06 

10  -0. 10000E+03 


0.0 

0.0 

0.  3  7500E»02 
0.37500F*02 
•0. 3750QF*02 
-0.375 00F*02 
0. 100Q0E*Q3 
0.100006*03 
-0. 100006*03 
-0.10000E*03 


0.200006*03 
Q.20000E*03 
0.1 0000E*03 
0.100006*03 
0.100006*03 
0.1 3000E*03 
0.0 
0.0 
0.0 
0.0 


COORDINATE  OESIGN  VARIABLES 
OESIGN  VARIABLE 


JOINT 


4 

5 

6 
7 
9 
9 

10 


X 

0 

0 

l 

l 

1 

1 

4 

4 

4 

4 


0 

0 

2 

2 

2 

2 

5 

5 

5 

5 


0 

0 

3 

3 

3 

3 

0 

0 

0 

0 


X 
0.0 
0.0 

-0. 10C0E*0i 
0.  10COE*01 


MULTIPLIER 

Y 

0.0 

8ll000E*01 
0.10006*01 


0.10005*01  -0.10005*01 
-0. IOOOF*01  -0. 10006*01 
-0.10006*01  0.10006*01 
0. 10006*01  0.1  GOOE  *Ci 

0.10006*01  -0.10006*01 
-0. 10006*01  -0.10006*01 


l 

0.0 

0.0 

0.10006*01 
0. lOCOE+Ol 
0.10006*01 
0.10006*01 
0.0 
0.3 
0.  0 
0.0 


MATERIAL  PROPERTIES 
r.QOE  6  RHO 


1  0.10006*08 


SIG-MIN 
0.1000E*00  -0.4000E+05 


S  IG-MAX 
0.  40006*05 


X. -EULER 
0 . 3927  E *02 
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Table  II 


2 5 -BAR  TOWER  SAMPLE  OUTPUT 

JOINT  .DISPLACEMENTS  CALCULATED  BY  SYSTEM  routine  "SAD" 


load  condition 

DISPLACEMENTS 
joint  „  l 

L  -0.219056- 

2  0.2 19106- 

3  3.907896- 

4  0.912786- 

5  -0.907896- 

6  -0.912786- 

7  0.0 

8  0.0 

9  0.0 

10  0.0 


JR  EE  06  FREEDOM 

0.3801 7E+00  - 

-0.38017E  +  00  - 

-0 . 1 5964E-0 1  ' 

0.17510E-0L 
0.159646-01  * 

-0.  1751 1E-0L 
0.0 
0.0 
0.0 
0.0 


■0.27099 E-OL 
■0.270996-01 
■0.687526-01 
0.36L  00E-01 
■0. 687526-01 
0.36100  6-01 
0.0 
0.0 
0.0 
0.0 


LOAD  CONDITION 
DISPLACEMENTS 


JOINT 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 


0.20L26E- 
0.229116- 
0.995316- 
J. 647336- 
0.814936- 
0.66536E- 
3.0 
0.0 
0.0 
0.0 


DEGREE  OF  FREEDOM 
2 

31  0 .3  88606  TOO  - 

)l  0.388606*00  - 

03  0.259516-01  ■ 

51  3:i«?Jfc8l  ' 

32  0.25192E-01 

0.0 

0.0 

0.0 

0.0 


C. 210 23 E- 01 
■3  .326876-0  L 
0.956536-01 
0  •  l Q297E+0Q 
0.62374E-01 
0.701946-01 
0.0 
0.0 
0.0 
0.0 


JOINT  COORDINATES 
JOINT  X 

1  -0.375006*02 

2  0.375006*02 

3  -0.375006*02 

4  0.375006*02 

5  0.375006*02 

6  -0.375006*02 

t  -0.100006*03 

3  0.100006*03 

9  0.100006*03 

10  -0,100006*03 

JOINT  COORDINATES 
JOINT  X 

L  -0.375006*02 

2  0.375006*02 

3  -0.375006*02 

4  0.375006*02 

5  0.375006*02 

6  -0.375006*02 

7  -0. 10000E*03 

8  0.100006*03 

9  0.100006*03 

10  -0.100006*03 


0.0 
0.  0 

0.375006*02 
0. 3750CE*02 
-0.375  006*02 
-0.375006*02 
0.  100005*03 
0.100006*03 
-0.  100006*03 
-0. 100 006*03 


0.0 
0.  0 

0.3  75  00E  +  02 
0.375006+02 
•0.375  006  +  02 
•0.375  006  +02 
0.  100  006+03 
0.  i 00006*03 
-0.100006*03 
-  0.  100  006*03 


0.2000UE+03 
0. 200006*03 
0.100006*03 
0.10000E+03 
0. 100006+03 
0. 100006*03 
0.0 
0.0 
0.0 
0.0 


0.200006+03 
0.2  0000E+C3 
0.100006+03 
0.10000E+03 
0.100006*03 
0. 100006+03 
0.0 
0.0 
0.0 
0.0 
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Table  III 


2 5 -BAR  TOWER  SAMPLE  OUTPUT 


AREA/COORDINATE  VARIABLES 


THERE  ARE  8 

INITIAL  VALUES 
0.200006 >01 


LOWER  HOUNDS 
0.10000e-04 


UPPER  BOUNOS 
0. 500006*02 


there  are  5 

INITIAL  VALUES 
0.375006*02 


LOWER  BOUNOS 

o.ioooce*oo 


UPPER  BOUNOS 
0.500006*03 


AREA  OESISN  VAR IABLES 

0.200006*01  0 .2  OOOOE*  01  0.200006*01  0.20000E+01  0.200006*01 

0.10030E-04  0.10000E-04  0.10000E-04  0.10000E-04  0.10000E-04 

0.50000E*02  0 .50000E*02  0.50000E+02  0.50000E+02  0.530006*02 
COORDINATE  VARIABLES 

0.375006*02  0.100006*03  0.l0000E*O3  0.10000E*03 

0.10000E*00  0.10000E+00  0.100006*00  0.1 OOOOE *00 
0.500006*03  0.50000E*03  0.50000E+03  0.500006*03 


JOINT  DISPLACEMENT  CONSTRAINTS 


OIPECT ION  i*X,  2*Y,  3*2,  0*RESULTANT 


LOAD 

LOWER 

UPPER 

NOOE 

OIR. 

CONO. 

BOUND 

BOUND 

l 

1 

1 

-0.35306*00 

0 • 35  0CE*O0 

2 

1 

l 

-0.35006*00 

0. 35  0CE*00 

1 

2 

1 

-0.35006*00 

0.35006*00 

2 

2 

1 

-0.35006*00 

0.3500E*00 

FREQUENCY  CONSTRAINTS 
LOWE"  BOUNOS  IN  CPS 

FREQUENCY  LOWER 

NUM86R  BOUND 

1  0.16006*02 
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Table  IV 


2 5 -BAR  TOWER  SAMPLE  OUTPUT 


FINAL  OPTIMIZATION  INFORMATION 


FINAL  optimization  information 


03J 


0.267256c *0  3 


DECISION  VAR148L6S  (X-VECTOR) 

1)  0.145956*91  0.74285E*00  0.983926*00 

7)  0.859S3F*)Q  0.69?  716*00  0.393676*02 

13)  0.682806*02 


CONSTRA 
l) 

71 
13) 

19) 

25) 

311 
37) 

43) 

49) 

55) 

61  ) 

67) 

73) 

79) 

85  1 
91  I 
971 
103) 
109) 
115) 

121  I 
127) 
133) 
139) 

1*5 ) 
151) 
157) 


INT  VALUES  <G 
”0.2  3806?*  91 
-0.292526-01 
-0.10090E+01 
—0.896076*  00 
-3.91922E* )0 
-0.102596*01 
-0.104916*01 
—0.754476  +00 
~0.1 102 lc  ♦  01 
-9.766  /OF* JO 
-0.111456*01 
-0.103666*01 
-0.195896*91 
-0. I 10185*01 
-0.923646+00 
-0.117266+01 
-0.806/76*00 
-0.115186*01 
-0 • 7  86026*03 
— C.84000E*00 
-0.8H20E*  JO 
-0.11316E*01 
— C«11028E*0l 
—0.122  34F*0 l 
-0.  72096E*  00 
-0.6806 7E* JO 
-0. I 183 1E*01 


VECTOR) 
0.993176*00 

‘  ‘  "526-1 “ 


0.292526-01 
0. 99 1 036  +  JO 
0. 110396*01 
0.1 J8076*01 
0.9  74136*00 
0.950886*00 
0.124566*01 
0.397896*00 
0.  123  3 IE* 01 
0. 385466*00 
0.  963  436*00 
0.941  HE  +  00 
0.8  98 1 7E  *00 
0.107646*01 
0. 327406*00 
0.  119326*01 
0.348156*00 
0.121406*01 
0.116006*01 
0.118886*01 
0.8  68  4  06 -*00 
0.897216*00 
0.776616*00 
0. 127906*01 
0.131 936*01 
0.316906*00 


THERE  ARE  7  ACTIVE  CONSTRAINTS 
CONSTRAINT  NUMBERS  ARE 

1  7  8  111  120  126  153 


THERE  ARE 


0  VIOLATED  CONSTRAINTS 


THERE  ARE 


0  ACTIVE  SIDE  CONSTRAINTS 


TERMINATION  CRITERION 
I TFR  EQUALS  ITMAX 


NUMBER  OF  ITERATIONS  »  20 

OBJECTIVE  FUNCTION  HAS  EVALUATED 
CONSTRAINT  FUNCTIONS  HERE  EVALUATED 
GRADIENT  OF  OBJECTIVE  HAS  CALCULATED 
GRADIENTS  OF  CONSTRAINTS  HERE  CALCULATED 


0. 1 J425E+01 
0.46934E+02 


-0.100686*01 
-0. 19707E+01 
-0.100356+01 
-0.63131E+00 
-0. 86905E ♦ 00 
-O.10420E+01 
-0. 10  797  E  +  01 
-0.849 20E *00 
— 0.106 27E+01 
-0.  856 83  E  +  00 
-0. 107  036*01 

-Ol ioioiltoi 

-0. 1 30666 ♦ 01 
-0.  969  89  F.-Ol 
-0.  1  709oc*  01 
-0.1 7345 E-04 
-0.38325E+00 
—0.272 19E*00 
-0. 15073  E*Ol 
-0.  1 3962  E*  01 
-0. L6847E+01 
-0.144736*00 
-0.21237E-01 
-0. 156126*01 


60 

60 

20 

20 


TIMES 
TIMES 
TIMES 
TI  ME  9 


0.11271E+01 
0. 13  8956  +03 


-0.1006dE*01 
*0.10 160E  +01 
■0.71 to5E *00 
-0.123496*01 
-0. 12349E+01 
-0.714656*00 
-0.121136*01 
-0. o55 JoE  *C0 
-0.O5536E+00 
-0.  12 11  36  +  01 
-0.104786*01 
-0.10478E+01 
-0. 101706*01 
-O.lOl/OE+Ol 
-0.110106*01 
-0.85  724 E  +00 
-0. 857246*00 
•0. 11 010E+01 
•J. 11 95  IE  +01 
-0.748776*00 
-0. 748776*00 
-0.11951E+01 
-0.994166*00 
■  0.909206  *00 
-0.99416E+00 
-0. 909196*00 


0.73567E*00 
0.  73636E  *02 


-0. 993166*00 
-0.98  405  E  *00 
-0. 12854E+01 
-0.765076*00 
-0. 765076*00 
-0.123546+01 
-0.78  8  70E  +00 
-0.  I  3*4oc ♦ 01 
-0.134466*01 
-0. 78870E+00 
-0.952166*00 
-0.952166*00 
-0.983046*00 
-0.983046*00 
-0. 399026*00 
-0.114286+01 
-0.114286+01 
-0.399026*00 
-0.804936*00 
-0.125126*01 
-0. 1 25 12E+01 
-0.80493E+00 
-0.10053E+01 
-0.109086+01 
-0.100586+01 
-0.109086+01 


•0. 197076*01 
-0.100636*01 
-0. 536846*00 
-0. 1 3813E+01 
-0.138136*01 
-  0. 5368*c*00 
-0.112986*01 
-0. 7oB32E*00 
-0. 788326*00 
-0.1 12986*01 
■  0. 1 04 1 26  +0 1 
*0.104 126 +01 
-0.10095E+01 
-0.10095E+01 
-0.147196+01 
-0. 332896*00 
-0. 33284E*00 
-0. 14719E+J1 
■0.175206*01 
-0.315506-01 
-0.315516-01 
•0.  17520E+01 
-0.982096*00 
■0. 721686*00 
-0. 98209E+0Q 
-0. 721586*00 
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V.  USER  GUIDE 


A.  INTRODUCTION 

In  developing  any  computer  code  for  engineering  analysis, 
it  is  necessary  to  additionally  develop  concise,  easily  under¬ 
stood  software.  This  USER  GUIDE  is  writen  to  be  easily 
followed  assuming  minimal  FORTRAN  knowledge.  The  format 
follows  that  of  the  optimization  code,  COPES/CONMIN , 
reference  1. 

This  chapter  is  devoted  to  acquainting  the  user  with  the 
code  and  necessary  input  data.  A  simple  3-bar  truss  analysis 
is  used  as  the  example . 

B.  GENERAL  FORMULATION 

Each  data  card  is  set  up  to  show  the  input  data  necessary 
with  the  3-bar  truss  example  underneath.  Formats  are  of  "15" 
and  "F10.0"  type.  "I"  formats  must  be  right  justified,  and 
"F"  formats  must  have  the  decimal  point.  The  number  of  cards 
read  per  data  block  is  listed  at  the  bottom  of  each  block. 

C.  CONSTRAINTS 

Constraints  are  calculated  and  stored  in  the  G  vector  as 
listed  in  the  following  chart.  The  total  number  of  constraints 
NCON  =  2 *NDSPLC  +  NFREQ  +  NE*NLC.  When  any  of  the  constraints 
are  missing  from  the  G  vector,  all  constraints  are  moved  up. 
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For  example,  if  there  is  no  frequency  constraint,  then  a  dis¬ 
placement  constraint  would  fill  the  first  location  of  the  G 
vector . 


r*. 


NFREQ 

(FREQUENCY  CONSTRAINTS) 


^1 

G1 


1 


1  2*  NDSPLC  *  NLC 

(DISPLACEMENT  CONSTRAINTS) 


all 


a 


1 


a 


NE 


a 


r 


NE 


1 

I 


MEMBER  #1 
LC  #  1 


1 


1 


MEMBER  #2 
LC  #1 


'J 


MEMBER  #NE 
LC  ¥NLC 


STRESS  AND 
BUCKLING 


STRESS  AND 
BUCKLING 

STRESS  AND 
BUCKLING 
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D.  EXAMPLE 


The  initial  layout  of  the  3-bar  truss  is  shown  in  Figure 
6.1.  Stress  constraints  were  imposed  as  well  as  constraints 
on  Euler  buckling,  displacement,  and  first  fundamental 
frequency . 

With  geometry  specified  as  per  the  figure,  2  independent 
load  conditions  (P1,P2)  were  imposed  with  3  member  sizing 
variables  (Al,  A2 ,  A3)  linked  so  that  A1=X(1),  A2=X(2), 

A3=X(1) . 

1 .  Properties /Conditions 

.  3 

The  material  used  had  a  density  (p)=0.1  lb/m  . 

7 

Young's  Modulus  was  selected  as  10  psi. 

LOADS:  P1=P2=  20,000  lbs. 

STRESS  LIMITS:  -15,000  <  o.  <  20,000  psi  i=l,3 
Symmetry  weis  maintained  by  linking  of  variables  so 
A1=A3.  Initial  design  began  with  A1=A2=A3=  1.0  in.  One  non- 
structural  mass  of  500  lbs  was  attached  at  joint  2.  Displace¬ 
ment  constraints  were  imposed  between  -.02  and  .02  inches  in 
the  x-direction  with  the  first  natural  frequency  limited  to 

a  value  w  <  20  Hz . 
n 

2  .  Input  Control  Parameters 

The  following  input  control  parameters  are  given  for 
ease  of  the  following  example.  NE=3,  NJ=4,  NCJ=4,  NMT=1,  NLC=2, 
MEIG  =  2 ,  NEIG1  =  2 ,  NFMASS  =  1,  NEJBC  =  1,  NDSPLC  =  2,  NFREQ  =  1 ,  LMASS  =  0 , 
IDVCLC=2.  Table  V  is  a  listing  of  commonly  used  nomenclature. 
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TiBLE  V 


COB  HON  VARIABLE  NOMENCLATURE 


A-  MEMBER'S  CROSS-SECTIONAL  AREA 
BL-  LOWER  BOUND  ON  DISPLACEMENTS 

BO-  OPPER  BOUND  ON  DISPLACEMENTS 

DIR-  DIRECTION  1=X,  2=Y.  3*Z, 

0*LIMIT  ON  RESULTANT  DISPL. 

E-  YOUNGS  MODULUS 
ELNO-  ELEMENT  NUMBER 

EPSEIG-  CONVERGENCE  TOLERANCE  OP  EIGENVALUE 
SOLUTIONS  (DEFAULT*. 0001) 

FC 1 ,  FC2  ,.  ...  PC  N-  LOWER  BOUND  ON  FIRST,  SECOND,  ETC. 

SYSTEM  NATURAL  FREQUENCY 
FX-  LOAD  PORCES  BEING  APPLIED  IN  THE  X  DIRECTION 

FY-  LOAD  FORCES  BEING  APPLIED  IN  THE  Y  DIRECTION 

FZ-  LOAD  FORCES  BEING  APPLIED  IN  THE  Z  DIRECTION 

GRAV-  ACCELERATION  DUE  TO  GRAVITY  (DEFAULT*  386. U) 
IX-  CONSTRAINT  IDENTIFIER.  IF  NON-ZERO  THE  X-DOF 
IS  CONSTRAINED 

II-  CONSTRAINT  IDENTIFIER.  IF  NON-ZERO  THE  Y-DOF 
IS  CONSTRAINED 

IZ-  CONSTRAINT  IDENTIFIER.  IF  NON-ZERO  THE  Z-DOF 
IS  CONSTRAINED 
JN-  JOINT  NUMBER  (GLOBAL) 

KEULER-  EULER  BUCKLING  COEFFICIENT  FOR  BAR 
ELEBENTS 

LC-  LOAD  CONDITION 

LM  ASS-  LUMPED  MASS  OPTIONS.  IF  LMASS  .  NE.  0  LUMPED 
MASS  MATRIX  USED 

IF  LMASS  *0  CONSISTENT  MASS  MATRIX  USED 
NEIG  1-NUMBER  OF  EIGENVALUE /EIGENVECTOR  TO  BE 
EVALUATED  DEFAULT  =  MIN(2*NEIG,  NEIG  +  8) 

NEUBC-  BUCKLING  CONSTRAINT  IDENTIFIERS  .  IF  NEUBC. 
NE.  0  EULER  BUCKLING  WILL  BE  IMPOSED  IN 
BAR  ELEMENTS 

NF MASS- NUMBER  OF  FIXED  MASSES  ATTACHED  TO  THE 
STRUCTURE 

NFREQ-NUHBEB  OF  FREQUENCY  CONSTRAINTS 

NID-  NUMBER  OF  INDEPENDENT  DEGREES  OF  FREEDOM 

NJ-NUMBER  OF  JOINTS 

NLC-  NUMBER  OF  LOADING  CONDTIONS 

NLJ-  NUMBER  OF  LOADED  JOINTS  FOR  THIS  LOAD  CONDITION 

NMT-  NUMBER  OF  SEPARATE  MATERIAL  TYPES 

RHO-  MATERIAL  DENSITY 

SIGHIN-  MINIMUM  ALLOWABLE  STRESS 

SIGMAX-  MAXIMUM  ALLOWABLE  STRESS 

XA-  INITIAL  VALUE  OF  AREA  DESIGN  VARIABLE 

XAL-  LOWER  BOUND 

XA U-  UPPER  BOUND 

XC-  INITIAL  VALUE  OF  COORDINATE  DESIGN  VARIABLE 
XCL-  LOWER  BOUNDS 
XCU-  UPPER  BOUNDS 
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The  following  user's  manual  is  divided  into  blocks  A 
through  M.  Appearing  directly  below  each  data  field  line 
are  the  parameters  for  the  3-bar  truss  example .  It  is  impor¬ 
tant  to  note  the  user  may  choose  any  units;  however,  all  units 
must  remain  consistent  throughout  the  problem. 
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£111  Si2£JS  i 

DESCRIPTION:  Title  Card 

Format  and  Example 


FIELD  CONTESTS 

1  AH  T  80  CHARACTER  TITLE 

MAT  BE  GIVEN  ON  THIS  LINE 
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slock  a 

DE  SCRIPT  ION:  Control  Parameters 

Format  and  Example 


NE 

NJ 

NCJ 

N»T 

NLC 

_ 

NEUBC 

_ 

NDSPLC 

NPRI 

_ 

FORMAT 

815 

pm 

■BHI 

■■■■ 

mmmmm 

H 

m 

i 

2 

...... 

i 

L 

2 

■  

0 

_ 

NEIG1 

NFHASS 

NFREQ 

IDTCLC 

IRECP 

FORMAT 

n 

2 

,  -  ..... 

1 

1 

2 

0 

L 

515 

i 

i. - — - — i,. . . ..  .t.  . . .  i 


NOTE:  DEFINITIONS  OF  INPUT  CONTROL  PARAMETERS  FOR 
PROGRAM  NEXT  PAGE. 
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ZIZJtR  SQSIiSI  | 

1  NE-number  of  elements 

2  NJ-nu*ber  of  joints 

3  HCJ-number  of  constrained  joints 

4  NMT-number  of  seperate  material  types 

5  HLC-number  of  load  conditions 

6  HEiJBC- buckling  constraint  identifier 

(If  N2UBC.NE.0  -EULER  buckling  con¬ 
straints  will  be  imposed  on  bar  ele.) 

7  NDSPLC -number  of  displacement  constraints 

8  NPRl-input  print  control  (if  NPBl.ne.O 

input  information  will  not  be  printed) 

1  NEIG-nuaber  of  precise  eigenvalues  to 

be  evaluated 

2  BEIG 1-number  of  eigenvalues  to  be  evaluated 

DEPAOLT=min.  Of  (2*NBIG  ,  NNEIG  +  8) 

3  HFHASS -number  of  fixed  masses  attached  to 

structure 

4  HFREQ- number  of  frequency  constraints 

5  iDVCLC-design  variable  control  parameter 

If  (I DVCLC.EQ.  1)  NDV*NDVAR 1 

If  (I DVCLC.EQ. 2)  NDV-KDVAR 1 ♦  NDVAR2 

If  (I DVCLC.EQ. 3)  HDV=HDVAR2 

6  IRECP-recripocal  variable  identifier 

If  (IR ECP.GT.O)  the  X-vector  contains 
the  rscripocal  of  area  and  VLB  £  VDB 
are  changed  accordingly 
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DA  21  JgiQCK  C 

DESCRIPTION;  Dynamic  Analysis  Information 


Format  and  Example 


LHASS 

_ 

GRAY 

EFSEIG 

FORMAT 

15, 

2F10.0 

0 

386. « 

0. 

. 

mu 

i 


2 


sauais 

LHASS -lumpe d  mass  option  (if  LHASS. NE.O) 
the  lumped  mass  matrix  is  used. 


GBAY-accleration  due  to  gravity 
(defau^aSSS  .'4  mches/sec  ) 


3 


EPSBIG-conv erge nee  tolerance  on  eigenvalu 
solution,  (default*. 0001) 


£114  block  d 

DESCRIPTION:  Joint  Coordinates  5 

oesiss.  laEiafele  Lialiaa 

Format  and  Example 


Z  jixjll  IZ  PCX  PCX 


PCZ  FORMAT 
15, 

3F1 0 . 0, 

315, 

3F10.0 


1  -10. 


0  1. 


FIELD 

CONTENTS 

1 

JN- joint  coordinate  number 

2 

X-x  coordinate 

3 

Y-y  coordinate 

U 

Z-2  coordinate 

5 

IX-design  variable  associated 

with  x  coord. 

6 

II-design  variable  associated 

with  y  coord. 

7 

IZ-design  variable  associated 

with  z  coord. 

8 

PCX -participation  coefficient 

of  x-coord. 

9 

PCI-partici pati  on  coefficient 

of  y-coord. 

10 

PCZ-partici pati on  coefficient 

of  z-coord. 

NOTE: 

L_ 

Nuaber  of  cards  read=NJ 
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UU  fiiO£K  E 


DESCBIPTION:  Material  Properties 

Format  and  Example 


E  ] 

_ 1 

BHO  j 

SIGHIN  ] 

_ ! 

SIGMAX 

KEOLER  1  FOBMAT 

J  5F10.0 

1000000.  | 
_ _ _ 1 

•  1 

-15000. 

20000. | 
t 

4. 

t  ■  1 

FIELD  COSTS  NTS 


1  E-Ioung*  s  Modulus 

2  BHO- material  density 

3  SIGM IN-ainiaum  allowable  stress 

4  SIGN AX-maximum  allowable  stress 

5  KEOLER-Euler  buckling  coefficient 

for  bar  element 


NOTE:  Number  of  cards  read»NHT 


PITA  BLOCK  F 

DESCRIPTION :  Bar  Element  Information 


Format  and  Example 


ELNO  NODE  1  ]  NODE  2  MAT  COD|  NDSG  |  A  | FORMAT  I 

!  i 515, F 1 0 . 0 ! 


I  1 


12  1 


3  2  1 


FIELD 


CONTENTS 


ELNO-element  number 
NODEl-global  number  associated  with 
element  node 

NODE2-global  number  associated  with 
element  node 

MATCOD -mate rial  type  of  this  element 
NDSG-design  variable  number  associated) 
with  this  element 
A-member  cross-sectional  area 


NOTE:  Number  of  cards  read*ME 
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DAJA  BLOCK  H 

OBIT  this  card  if  NL  C=0  vas  read  in  block  B. 
DESCRIPTION;  Joint  Loading  Information 

Format  and  Example 


r  nlj  r 

i  i 

i _ i 


i — r 


FORMAT  I 


1 


15 


IZ 


JN  I  FX 

j 

FT 

FZ 

_ 

[format 

1 15 , 3F  1 0 . 0 

.1 _ 


r1 

14140. 

_ _ 

-14140. 

_ 

ZH 

, _ 

- 1 

1 

_ 1 

-14  140. 

-14140. 

0. 

. . 

L  _ — 1 

1 

L-  ,  .  1 

FIELD 


CO  NT  a?  T 


1  NLJ -number  of  loaded  joints 

for  this  load  condition 

1  JN- joint  number 

2  FX 

3  FT-  Forces  in  the  X,Y,Z  directions 

4  FZ 
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Dili  BLOCK  I 

OBIT  this  block  if  NFMASS  =  Q  was  read  in  block  B. 
DESCRIPTION;  Lumped  Mass  Information 

Format  and  Exaaple 


JN 

_ - _ 

HASS 

_ 

FORMAT 

I5,F10 . 0 

LLZ 

500. 

— 

_ 

FIELD  CONTENT? 

1 
2 


JN-joint  number 

MASS -concentrated  mass  at  joint  (JN) 
mass  is  in  force  units 


I 


I 

f 


i 


1 
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DA2 A  BLOCK  J 

OMIT  this  block  if  NDVAH1=0 

DESCSIPTION;  Design  Variable  Information 

(AREA  Variables) 

Format  and  Example 


FIELD 


CONTENTS 


1  XA- initial  value  of  area  design  variables 

2  XAL -lower  bounds  on  area  design  variables 

3  XACJ-upper  bounds  on  area  design  variables 

NOTE:  read  one  value  of  XA,  XAL,  XAO  for  each 

independent  area  variable  defined  in  Block  D 


Number  of  cards  read  =as  required 
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DATA  BiO£K  K 

OMIT  this  block  if  NDVAR2  =  0 

DESCRIPTION:  Design  Variable  Information 

(COORDINATE  Variables) 


Format  and  Example 


— 

XC  (1) 

xc<2)  :  ...  . 

XC  (NDVAR1) 

FORMAT 

"1 

8F10.0 

_ 

- 1 — - 

i  10.  | 

L _ 1 

- r  “i - 1 

_ _ L  J 

r  i.o  ] 
1 _ l 

i  1  1 XCL 

_ I  _J  J 

r 

1  ieD 

£I2LM 

1 
2 
3 

NOTE:  read  one  value  of  XC,  XCL,  XCD  for 

each  independent  coordinate  variable 
defined  in  Block  D 


CONTENT? 

XC-initial  value  of  coord,  design  variables 
XCL- lower  bounds  on  coord,  design  variables 
XCU-upper  bounds  on  coord,  design  variables 


Number  of  cards  read  -as  required. 


DATA  BLOCK  L 

OMIT  this  block  if  NDSPLC*0  in  Block  A 
DESCRIPTION;  Joint  Displacement  constraint  Information 

Format  and  Example 


DIR 

LC 

BL  | 

BO  j 

FORMAT  j 

f 

1 

L _ J 

1 

I 

I 

3I5,2F10.0 
_ 1 

1  2  i 

J  J 

rrn 

1 

-.02 

.02  ] 

j 

L _ 1 

2 

1  1 
_ 1 

-.02] 
L _ J 

I  .02 

L 

1LSU2 


CONTJNIS 


1  JN-joint  number 

2  DIR-i irection  1  *X  2=T  3=z 

0=liait  on  resultant  displ. 

3  LC-load  condition 

4  BL-lower  bound  on  displacement 

(If  DIR*0  read  0  here) 

5  BO-upper  bound  on  displacement 


Number  of  cards  read*  NDSPLC 


DAJA  BLOCK  H 

OfllTthis  block  if  NFREQ=0  was  read  in  Block  A 
DESCRIPTION;  Frequency  Constraint  Inforaation 

Foraat  and  Example 


FIELD  CONT  ENTS 

1  FC1-  lower  bound  on  first  natural  frequency 

constraint  in  CPS.  (cycles  per  second) 
N  FCN-  lower  bound  on  NFREQ-th  natural 

frequency  constraint  in  CPS. 

(cycles  per  second) 


NOTE;  OHIT  this  block  if  NFREQ*0  was  read  in  data 
block  B. 

Number  of  cards  read3  as  required 
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VI.  NUMERICAL  EXAMPLES 


A.  INTRODUCTION 

Design  of  planar  trusses  and  3-dimensional  space  towers 
are  presented  here  and  the  corresponding  numerical  results 
are  summarized  to  demonstrate  the  use  of  the  code. 

In  the  examples  given  here,  the  design  variables  are 
member  cross-sectional  areas  and  joint  coordinates.  It 
should  be  emphasized  that  in  practical  design,  the  reciprocal 
of  the  member  areas  is  usually  a  better  choice  for  the  design 
variables.  However,  the  purpose  here  is  to  knowingly  create 
difficult  optimization  problems,  thus  the  choice  of  variables. 

The  examples  begin  with  the  3-bar  truss . 

B.  CASE  1:  3-BAR  PLANAR  TRUSS 

The  simple  3-bar  planar  truss,  as  shown  in  Figure  5.1  has 
been  previously  used  for  the  user  guide  example.  This  struc¬ 
ture  was  designed  for  optimum  geometry  subject  to  a  set  of 
two  load  conditions,  buckling  constraints,  displacement  con¬ 
straints,  and  a  lower  bound  on  the  system's  first  natural 
frequency.  The  allowable  stresses  specified  are 

-15000.  <  <  20000.  psi  i=l,3 

The  member  areas  were  linked  in  the  following  groups:  A1=A3; 
A2  . 
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1.  Case  la 


Only  the  stress  constraint  was  imposed  for  this  case. 
Final  design  information  is  given  in  Table  VI.  The  number  of 
analysis  for  this  design  was  38. 

2  .  Case  lb 

This  case  includes  the  previous  constraint  plus  con¬ 
straints  on  buckling  < KEULER  =  4 . o )  ,  displacement  (-.02  -.02  in 
and  first  natural  frequency  ((^=20.0) .  Results  are  given  in 
Table  VII.  The  number  of  analysis  required  for  this  design 
was  16  with  1.73  seconds  of  CPU  time. 


Figure  6.1  3-3AR  TRUSS 
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C.  CASE  2:  18-BAR  PLANAR  TRUSS 


A  cantiliver  truss,  as  shown  in  Figure  6.2,  has  been  pre¬ 
viously  used  as  a  standard  test  case  for  structural  design 
[Ref.  2].  The  structure  was  analyzed  for  a  single  set  of  load 
conditions  with  allowable  stresses  being 

-20000.  <  a ^  <  20000.  psi  i  =  l,13 

Young's  modulus  was  taken  as  10 ^  psi  with  a  material  density 
P=0 . 1  lb . /cu  in . 


Figure  6.2  18-BAR  TRUSS 
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The  independent  coordinate  variables  were  taken  as  X3,  Y3, 
X5,  Y5 ,  X7,  Y7,  X9 ,  Y9 .  The  member  areas  were  linked  as 
follows:  A1=A4=A8=A12=A16 ;  A2 =A6 =A10 =A14 = A1 3 ;  A3=A7=A11=A15 ; 

A5=A9=A13=A17 .  There  are  a  total  of  four  independent  area 
variables  and  eight  coordinate  variables. 

1 .  Case  2a 

This  design  analysis  includes  stress,  buckling  (KEULER= 
4.05,  displacement  (-10.0  to  10.0  in.),  and  a  first  fundamental 
frequency  of  3  Hz.  Additionally,  a  nonstructural  mass  of  W=50Q0 
lbs.  was  attached  to  node  1.  The  number  of  analyses  for  this 
design  was  62.  Results  are  presented  in  Table  VIII.  26.72 
seconds  of  CPU  time  were  used  for  this  calculation. 

D.  CASE  3:  2 5 -BAR  SPACE  TOWER 


Figure  6.3  25-BAR  SPACE  TOWER 
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The  25-bar  space  tower  shown  in  Figure  6.3  was  designed 
for  two  independent  load  conditions  given  in  Table  X.  The 
allowable  stresses  were  specified  as  follows 

-40000.  <  <  40000.  psi  i=l,25 

7 

Young's  modulus  was  selected  as  10  psi  with  a  material  density 
p=.l  lb/cu.  in.  Members  are  assumed  tubular  with  a  nominal 
diameter  to  thickness  ratio  of  D/t=100  giving  KEULER= 39 . 2 74 . 
Symmetry  was  imposed  in  both  the  x-z  and  y-z  planes.  Non- 
structural  masses  of  W=5G0  lbs.  were  attached  at  nodes  1  and 
2.  Coordinate  variables  were  X4 ,  Y4,  Z4,  X3,  and  Y3  with  the 
remaining  coordinates  linked  to  maintain  symmetry.  Area  vari¬ 
ables  were  linked  in  the  following  manner:  A1 ;  A2=A3=A4=A5; 
A6=A7=A8=A9 ;  A10=A11;  A12=A13;  A14=A15=A16=A17 ;  A18=A19 =A20 =A21 
A22=A23=A24=A25 . 

1 .  Case  3a 

Stress,  displacement  (-0.35  to  0.35  in.),  Euler  buck¬ 
ling  and  first  natural  frequency  limited  to  a  value  >  16 
Hz.  were  imposed  for  this  case.  Final  design  information  is 
given  in  Table  XI.  The  number  of  analyses  for  this  design  was 
61  with  27.99  seconds  of  CPU  time  used. 

E.  CASE  4;  234-BAR  SPACE  TOWER 

The  inital  layout  of  the  tower  is  shown  in  Figure  6.4 
stress  limits  were  as  follows 

-15000.  <  <  20000.  psi  i=l,234 

7 

Young's  modulus  was  chosen  as  steal  3.xl0  psi  with  a  material 
density  of  aluminum  p  =.l  lbs/cu.  in. 


1.  Case  4a 


Stress  constraints  as  well  as  constraints  on  Euler 
buckling  were  imposed.  The  resulting  problem  had  56  area 
design  variables  and  42  coordinate  variables.  This  problem 
required  91  structural  analyses  using  72  minutes  of  CPU  time. 
The  objective  function  was  evaluated  89  times  and  gradients 
were  calculated  30  times.  Results  are  presented  in  Tables  XII 
through  XIV.  Although  the  weight  of  the  structure  was  only 
reduced  by  2%,  it  should  be  noted  that  constraints  were  ini¬ 
tially  violated  and  the  optimizer  overcame  the  constraint 
violations.  It  is  believed  the  optimum  that  was  reached  is 
not  the  true  optimum  because  of  the  extreme  non-linearity  of 
the  problem. 
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left/right 


SPACE  TOWER 


Table  VI 


3- 

BAR  TRUSS 

.  DESIGN  INFORMATION 

(STRESS) 

BODY 

LOAD 

JOINT 

LOADS 

COND. 

X 

~^T~  Z 

1 

2 

0 . 1414E  +  0  5 

-0.1414E+05  0.0 

2 

2 

-0 .1414E+0S 

-0 . 1414E  +  0  5  0.0 

FIXED  MASSES 

NO. 

JOINT 

MASS 

1 

2 

500  lbs. 

MEMBER 

AREA  (sq.  in.) 

INITIAL 

FINAL 

A1=A3 

1.0 

.  7662E  +  00 

A2 

1.0 

.4833E+00 

JOINT 

COORDINATES  (in.) 

INITIAL 

FINAL 

2 

10 

.9  35  6E  +  0 1 

FINAL  WEIGHT 

FOR  AREA 

AND  COORDINATES 

WEIGHT  =  2.6349  (lbs.) 
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Table  VII 

3-BAR  TRUSS.  DESIGN  INFORMATION 
(STRESS,  DISPLACEMENT,  BUCKLING,  FREQUENCY) 

BODY 

JOINT  LOADING  INFORMATION 
SAME  AS  TABLE  VI 


JOINT  DISPLACEMENT  INFORMATION 


JOINT 

DIR. 

LOAD 

LOWER 

UPPER 

NUMBER 

COND. 

30UND 

BOUND 

2 

1 

1 

-.02 

.02 

2 

2 

1 

-.02 

.02 

AREA  (sq.  in. 

l 

MEMBER 

INITIAL 

FINAL 

A1=A3 

1.0 

0 .99995E+00 

A2 

1.0 

0.1000  E-01 

COORDINATES  (in.) 

JOINT 

INITIAL 

FINAL 

2 

10  .0 

0  . 999  86E  +  01 

FINAL  WEIGHT  =  2 . 9 3 81  ( lbs . ) 

CRITICAL  CONSTRAINT  -  FREQUENCY  j 

I 
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Table  VIII 


13-BAR  TRUSS. 

LOAD 

CONDITIONS 

BODY 

LOAD  CONDITIONS 

(lbs  . ) 

LOAD  COND. 

JOINT 

LOADS 

X 

Y 

Z 

1 

1 

0.0 

-0 .2000E+05 

0.0 

1 

2 

0.0 

-0 .2000E+05 

0.0 

1 

4 

0  .0 

-0 . 2000E+05 

0  .0 

1 

6 

0.0 

-C . 2000E  +  0  5 

0  .0 

1 

8 

0.0 

-0 .2000E+05 

0.0 

FIXED  MASS  INFORMATION 

NO. 

JOINT 

MASS 

1 

1 

5000  (lbs 
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Table  IX 

13-BAR  TRUSS.  DESIGN  INFORMATION 
(STRESS,  DISPLACEMENT,  EULER  BUCKLING,  FREQUENCY) 


AREA  (sq.  in.) 


MEMBER 

INITIAL 

FINAL 

A1=A4=A8=A12=A16 

10.0 

38.20 

A2=A6=A10=A14=A18 

10.0 

38.27 

A3=A7=A11=A1S 

10 .0 

15.98 

A4=A9=A13=A17 

10.0 

20.62 

COORDINATES  (in.) 

JOINT 

INITIAL 

FINAL 

X  Y  X  Y 


3 

1000  . 

0  . 

859  . 

0  . 

5 

750  . 

0  . 

628  . 

0  . 

7 

500  . 

0  . 

384  . 

0  . 

9 

2  50  . 

0  . 

148  . 

0  . 

FINAL  WEIGHT  =  13,940  (lbs.) 
CRITICAL  CONSTRAINT  -  FREQUENCY 
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Table  X 


25-BAR  TRUSS.  LOAD  CONDITIONS 

30DY 


LOAD 

CONDITIONS  (lbs.) 

LOAD  COND. 

JOINT 

X 

LOADS 

Y 

1 

1 

0.0 

2000  . 

2 

0.0 

-2000  . 

2 

1 

1000  . 

10000  . 

2 

0  . 

10000  . 

3 

500  . 

0  . 

6 

500  . 

0  . 

FIXED  MASS  INFORMATION 

NO.  JOINT  MASS 

2  1  500  . 

2  500  . 


Z 

-5000  . 
-5000  . 
-5000  . 
-5000  . 

0  . 

0  . 


(lbs  .) 
(lbs  .) 
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Table  XI 


25-BAR  TRUSS.  DESIGN  INFORMATION 
(STRESS,  DISPLACEMENT,  EULER  BUCKLING,  FREQUENCY) 


AREA  (sq.  in.) 

MEMBER 

INITIAL 

FINAL 

A1 

2.0 

1.46 

A2=A3=A4=A5 

2.0 

0.74 

A6=A7=A8=A9 

2  .0 

0.98 

A10  =A11 

2  .0 

1.04 

A12  =A13 

2.0 

1.12 

A14=A15=A16=A17 

2  .0 

0.74 

A18=A19=A20=A21 

2  .0 

0.86 

A22=A23=A24=A25 

2.0 

0.70 

COORDINATES  (in.) 

JOINT 

INITIAL 

X  Y  Z 

X 

FINAL 

Y 

z 

4 

37.5  37.5  100 

39.4 

46.9 

138.9 

8 

100.0  100.0  0.0 

73.6 

68  .  3 

0.0 

FINAL  WEIGHT  =  267  (lbs.) 
CRITICAL  CONSTRAINT  -  FREQUENCY 
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Table  XII 


234-BAR 

SPACE 

TOWER.  LOADING 

INFORMATION 

BODY 

LOAD 

CONDITIONS  (lbs. 

2. 

LOAD  COND. 

JOINT 

LOADS 

X 

Y 

z 

1 

49 

6000  . 

-20000  . 

0.0 

1 

50 

6000  . 

-20000  . 

0.0 

2 

49 

6000  . 

-20000  . 

0.0 

2 

50 

-6000  . 

-20000  . 

0.0 

3 

49 

6000  . 

-20000  . 

0.0 

3 

50 

3000  . 

-10000  . 

5000  . 

4 

49 

3000  . 

-10000 . 

-5000  . 

4 

50 

3000  . 

-10000 . 

5000  . 

5 

49 

-3000  . 

10000  . 

5000  . 

5 

50 

-3000  . 

10000  . 

-5000  . 

FIXED 

MASS  INFORMATION 

NO. 

JOINT 

MASS 

1 

49 

200  . (lbs.) 

2 

50 

200. (lbs.) 
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Table  XIII 


234-BAR  SPACE  TOWER.  DESIGN  INFORMATION 


AREA  Csq.  in.) 


MEMBER 

INITIAL 

FINAL 

A1,A2 ,A3,A4 

25.0 

25.38 

AS ,A6 ,A7,A8 

25.0 

25  .  35 

A9 , A10 , All , A12 

25.0 

25.28 

Ak3,A14 ,A15 ,A16 

25.0 

25.21 

A17,A18,A19 ,A20 

25  .0 

25.14 

A21 ,A2  2 ,A2  3 ,A24 

25.0 

25.07 

A2 5  , A26  , A2  7  , A2  8 

25.0 

25.00 

A29  ,A30  , A31  ,A32 

25.0 

24.95 

A33,A34,A35 ,A36 

25.0 

24.89 

A37 ,A38 ,A39 ,A40 

25.0 

24.35 

A41,A42,A43,A44 

25.0 

24.98 

A45 , A46 ,A47,A48 

25.0 

24 .87 

A49,A50,A51,A52 

25.0 

24.88 

A53,A54,A55,AF6 

25.0 

24.89 

A5  7 , A5  8 ,A5  9 ,A6  0 

25.0 

24 . 89 

A6 1 , A6  2 , A6  3 , A6  4 

25.0 

24.90 

A6  5 ,A6  6 ,A6  7 , A6  8 

25.0 

24.91 

A69 , A70 ,A71,A72 

25.0 

24.91 

A73,A74,A75,A76 

25.0 

24.91 

A77,A78,A79,A80 

25.0 

24.91 

A81 , A8  3 

25.0 

24.96 

A82 , A84 

25.0 

24.96 

A85 ,A87 

25.0 

24.96 

A86 ,A88 

25.0 

24.95 

A89 ,A90 , A9 1 ,A9  2 ,A9  3 , A94 ,A9  5 ,A96 

25.0 

24.68 

A9  7 ,A9  8 ,A99 ,A10Q ,A101,A102 ,A103,A104 

25.0 

24.69 

A150 ,A106 ,A10  7 ,A10  8 ,A109 ,A110 ,A111,A112 

25  .0 

24.71 

A113 ,A114 ,A115 ,A116 ,A117,A118 ,A119 ,A120 

25.0 

24.73 

A121 , A12  2 ,A12  3 ,A124 ,A12  5 ,A126 ,A127,A128 

25.0 

24 .74 

A129 ,A130 , A1 31 , Al 32 ,A133 ,A134 ,A135 ,A136 

25.0 

24.75 

A137 ,A138 ,A139 , A140 , A141 , A14 2 ,A143,A144 

25.0 

24.74 

A14  5 , A146 ,A14  7,A14  8,A14  9 ,A150 ,A151,A15  2 

25.0 

24.72 

A15  3 ,A154 ,A155 ,A15 6 , A15 7 , A15 8 , A15 9 ,A16  0 

25  .0 

24  .69 

A161,A16  2 , A16  3 ,A164 , A16  5 ,A166 ,A167,A168 

25  .0 

24 .65 

A169 ,A170 ,A171,A172,A173,A174 ,A175 ,A176 

25.0 

24 . 82 

A177,A178,A179 ,A130 

25.0 

24.96 

A181 ,A182 ,A19  3 , A134 

25.0 

24.96 

A135 ,A186 ,A137 ,A138 

25  .0 

24  .  97 

A189 , A190 ,A191 , A19  2 

25  .0 

24.86 
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Table  XIII  (cont'd) 


AREA  (sq.  in.) 


MEMBER 

INITIAL 

FINAL 

A19  3 ,A194 ,A19  5 ,A196 

25.0 

24.95 

A197 ,A198 ,A199 ,A200 

25.0 

24.97 

A201.A202 

25.0 

24 .91 

A20  3  ,A204 

25.0 

24.92 

A205,A206,A207,A208 

25  .0 

24 . 86 

A209 , A2 10 ,A211 ,A2 12 

25.0 

24 .87 

A213 ,A214 

25.0 

24.89 

A215 , A2 16 

25.0 

24.90 

A2 17 ,A2 18 

25.0 

24.91 

A219 ,A220 

25.0 

24.92 

A221,A222 

25.0 

24.92 

A22  3 ,A224 

25.0 

24.93 

A225 ,A226 

25.0 

24.93 

A22  7  ,A22  8 

25.0 

24.94 

A2  2  9  ,A2  30 

25.0 

24.94 

A2  31 , A2  32 

25.0 

24.94 

A2  3  3  , A2  34 

25.0 

24.93 

COORDINATES 

C  in . ) 

JOINT 

INITIAL 

FINAL 

X 

Y 

Z 

X 

Y 

Z 

1 

120 

0* 

120 

120  .5 

0 

121.4 

5 

111 

120 

111 

110.7 

120.4 

112.1 

9 

102 

240 

102 

102.2 

241.2 

103.2 

13 

93 

360 

93 

93.1 

362  .6 

93.7 

17 

94 

4  80 

84 

34.0 

483.8 

34.4 

21 

75 

600 

75 

74 .9 

604.5 

75  .2 

25 

66 

720 

66 

65.9 

724.2 

65.9 

29 

57 

840 

57 

56  .  8 

842  .5 

56  .  3 

33 

48 

960 

48 

47 . 8 

960 .4 

47.3 

Initial  weight  * 

84,524  lbs. 

Final  weight  =  83 

,142  lbs. 

•Final 

weight  (AREA  V ARB I ABLE S 

ONLY)  => 

28,900  lbs. 

*  This  clearly  indicates  optimum  not  reached  with  previous  case. 
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Table  XIII  (cont’d) 


COORDINATES  (in.) 


JOINT 

INITIAL 

FINAL 

X 

Y 

Z 

X 

Y 

Z 

37 

39 

10  80 

39 

3  8.8 

1077.4 

38  .  8 

41 

30 

0* 

30 

29  .9 

- 

29.9 

45 

30 

1248 

30 

29  .  8 

1212.0 

29.9 

49 

0* 

1243 

90 

- 

1203.5 

90  .0 

51 

15 

1224 

60 

15  .0 

1194.1 

58.9 

53 

15 

1248 

60 

15.0 

1214.7 

59.2 

*NOT  A  DESIGN  VARIABLE 
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VII .  CONCLUSIONS  AND  RECOMMENDATIONS 

A.  CONCLUSIONS 

The  finite  element  code  was  presented  coupled  to  an 
optimizer  for  truss  analysis  and  design.  Trusses  were 
designed  for  minimum  weight  with  multiple  load  conditions 
considered . 

The  displacement  method  for  static  analysis  and  the  sub¬ 
space  iteration  method  ior  eigenvalues  were  applied. 

Several  examples  were  considered.  In  every  case  the  code 
worked  as  an  analysis  tool,  and  significant  weight  reductions 
were  obtained  with  the  coupled  optimizer  CONMIN.  Run  times 
and  results  compared  to  the  same  test  cases  with  other  codes 
indiate  that  the  code  is  competitive  as  a  design  tool. 

B.  RECOMMENDATIONS 

The  following  recommendations  may  be  of  value  for  follow 
on  work: 

1.  The  code  should  be  modified  so  the  user  can  access 
stresses  directly. 

2.  The  code  should  be  extended  to  include  other  elements 
such  as  frames  and  plates . 

3.  An  out  of  core  equation  solver  should  be  added. 

4.  The  method  of  gradient  calculation  should  be  dependent 
on  specific  gradients  required  reference  3  and  reference 
4  . 
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5.  Gradients  of  frequency  constraints  would  benefit  from 
a  more  efficient  algorithm  reference  5. 

6 .  The  need  for  a  large  scale  public  structural  optimiza¬ 
tion  code  still  exists  . 
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APPENDIX  A 


DATA  FILES 


A.  INTRODUCTION 

This  appendix  contains  the  data  files  used  to  create  the 
test  cases  in  Chapter  VI.  Additionally  the  data  file  for  the 
user’s  guide  in  complete  form  is  presented. 
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Table  XIV 


DATA  FILE  3-BAR  TRUSS 


(STRESS) 


3-  BAR  TRUSS  (TEST  CASE  A-STRESS  ONLY J 


3 

0 

0 

1 

2 

3 

4 


4  4 

0  1 

386.4 

-10. 


0. 

0. 


10. 

10000000. 


3 

4 


3 

4 
1 
0 


.1 

2 

2 

2 

1 

0 

1 

1 


l 
0 

0.0 

0. 

-10. 

0. 

0. 


0. 

0. 

0. 

0. 


-15000.  20000. 

I  l  !• 


1 

0 

0 

1 


0 

0 

0 

0 


0 

0 

0 

0 


2 

1 


1. 

1. 


Sol 

10. 

10. 

1.0 

20. 


14140.  -14140. 


40.  -14140. 

00. 


’“I 


0. 


i*. 
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j  j 

Table  XV  i j 

DATA  FILE  3-BAR  TRUSS  !  | 

(STRESS,  DISPLACEMENT,  BUCKLING,  FREQUENCY)  j 


3- 


BAR 

3 

2 

0 

1 

2 

3 

4 


10000000. 


TRUSS 
4  4 

2  1 

386.4 

-10. 

0. 

0. 

10. 


(TEST  CASE 
1 
1 

0.0 

0. 

-10. 

0. 


D-STRESS.DISPL. ,3UCKL. ,FREQ.  I 
2  12  0 
2 


L 

2 

3 

1 

2 

3 

4 


1. 

.01 

10. 

10. 

1.0 

20. 


1 

3 

4 
1 
0 
1 
1 


.1 

2 

2 

2 

h 

1 


0. 

1 

1 

1 

1 

1 

1 

1 


-15000. 
1 
2 
1 


14140.  -14140. 


-14140. 

500. 

1. 

.01 

10. 


-14140. 


20. 


2 

2 


1 

2 


1 

1 


-.02 

-.02 


0. 

0. 

0. 

0. 

20000. 

1. 

1. 

1. 


1 

0 

0 

1 


0 

0 

0 

0 


0 

0 

0 

0 


-1. 

0 

0. 

1. 


.02 

.02 


1 


i 


' 
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Table  XVI 


DATA  FILE  18 -BAR  TRUSS 
(STRESS,  DISPLACEMENT,  BUCKLING,  FREQUENCY) 

1S-RAR  TRUSS  CASE  {TEST  CASE  STRESS/O ISP- /SUCKLING/ FREQ. i 


13 

2 

0 


U  11 
4  1 

386.4 


1 

l 

0. 


1 

2 


i 

3 

1250. 

1000. 

1000. 

250. 

250. 

1 

2 

0  1. 

1. 

4 

5 

750. 

750. 

250. 

0. 

3 

4 

0  1. 

1 

l. 

6 

7 

5  00. 
500. 

250. 

0. 

5 

6 

1. 

1.  ( 

8 

9 

250. 

250. 

250. 

0. 

7 

8 

1. 

1.  ! 

10 

0. 

250. 

i 

11 

10000000. 
1 
2 

3 

4 

5 

6 
7 
3 
9 

10 
11 
12 
13 

15 

16 

17 

18 
1 
2 

3 

4 

5 

6 
7 

a 

9 

10 


lk 


l 

1 

2 
2 
3 

3 

4 

4 

5 

5 

6 
6 
7 

7 

8 
8 
9 
9 
0 


2 

3 

3 

4 
4 

I 

6 

6 

7 

7 

a 

9 

10 

!2 

0 


-20000. 

1  10. 

2 

3 
1 

4 
2 

3 
1 

4 
2 

3 
1 

4 


20000. 


-20000 

-20000 

-20000 


30. 

1. 

1000. 

1000. 

too. 

1200. 

1 

4.0 


5000. 

30. 

1. 

1000. 

0. 

-50. 

200. 

2  1 


-18888 


30. 

1. 

1000. 

750. 

500. 

1000. 

-10. 


10. 

10. 

10. 

10. 

10. 

10. 

10. 

10. 

10. 

10. 

10. 

10. 

10. 

10. 

10. 

10. 

10. 


30. 

1. 

^000. 

-50. 

150. 

10. 


503. 

*88: 


0. 

-50. 

150. 


ft0.' 

500. 


150.' 


I 
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Table  XVII 
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Table  XVIII 


DATA  FILE  2 34 -BAR  SPACE  TOWER 


3*-BAR 

SPACE  TOWER  DESIGN 

(TEST  CASE 

STRESS. BUCK. . 

DISPL..FREQ. 

) 

234 

58  4 

1  5 

1  20 

0 

1 

3  2 

1  2 

0 

386.4 

0. 

l 

120. 

0. 

120. 

1 

0 

2 

1. 

0. 

1. 

2 

-120. 

0. 

120. 

l 

0 

2 

-1. 

0. 

1. 

3 

-120. 

0.0 

-120. 

l 

0 

2 

-1. 

0. 

-1. 

4 

120. 

0. 

-120. 

1 

0 

2 

1. 

0. 

-1. 

5 

111. 

120. 

111. 

3 

4 

5 

1. 

1. 

1. 

6 

-III. 

120. 

111. 

3 

4 

5 

-1. 

1. 

1. 

7 

-111. 

120. 

-111. 

3 

4 

5 

-1. 

1. 

-1. 

a 

111. 

120. 

-111. 

3 

4 

5 

1. 

1. 

-1. 

9 

102. 

240. 

102. 

6 

7 

8 

1. 

1. 

1. 

1? 

-102. 

240. 

102. 

6 

7 

a 

-1. 

1. 

1. 

-102. 

240. 

-102. 

6 

7 

8 

-1. 

1. 

-1. 

12 

102. 

240. 

-102. 

6 

7 

8 

1. 

1. 

-1. 

13 

93. 

363. 

93. 

9 

10 

11 

1. 

1. 

1. 

14 

-93. 

360. 

93. 

9 

10 

11 

-1. 

1. 

1. 

15 

-93. 

360. 

-93. 

9 

10 

11 

-1. 

1. 

-1. 

16 

93. 

360. 

-93. 

9 

10 

11 

1. 

1. 

-1. 

17 

84. 

480. 

84. 

12 

13 

14 

1. 

1. 

1. 

18 

-84. 

480. 

84. 

12 

13 

14 

-1. 

1. 

1. 

19 

-84. 

480. 

-84. 

12 

13 

14 

-1. 

1. 

-1. 

20 

84. 

480. 

-84. 

12 

13 

14 

1. 

1. 

-1. 

ii 

75. 

600. 

75. 

15 

16 

17 

1. 

1. 

1. 

-75. 

600. 

75. 

15 

16 

17 

-I. 

1. 

1. 

23 

-75. 

600. 

-75. 

15 

la 

17 

-1. 

1. 

-1. 

24 

75. 

600. 

-75. 

15 

16 

17 

1. 

1. 

-1. 

25 

66. 

720. 

66. 

18 

19 

20 

1. 

1. 

1. 

26 

—66. 

720. 

66. 

18 

19 

20 

-1. 

1  <* 

1. 

27 

-66. 

720. 

-66. 

18 

19 

20 

-1. 

1. 

-1. 

28 

66. 

720. 

-66. 

16 

\l 

20 

1. 

1. 

-1. 

29 

57. 

340. 

57. 

21 

23 

1. 

1. 

1. 

30 

-57. 

840. 

57. 

21 

22 

23 

-1. 

1. 

1. 

31 

-57. 

840. 

-57. 

21 

22 

23 

-1. 

1. 

-1. 

32 

57. 

840. 

-57. 

21 

22 

23 

1. 

1. 

-1. 

33 

48. 

960. 

48. 

24 

25 

26 

1. 

1. 

1. 

34 

—48. 

9SQ  . 

48. 

24 

25 

26 

-1. 

1. 

1. 

35 

-48. 

960. 

-4b. 

24 

25 

26 

-1. 

1. 

-1. 

36 

48. 

960. 

-48. 

24 

25 

26 

1. 

1. 

-1. 

37 

39. 

1080. 

39. 

27 

28 

29 

1. 

1. 

1. 

38 

-39. 

1380. 

39. 

27 

28 

29 

-1. 

1. 

1. 

39 

-39. 

1080. 

-39. 

27 

28 

29 

-1. 

1. 

-1. 

40 

39. 

1030. 

-39. 

27 

28 

29 

1. 

1. 

-1. 

41 

30. 

1200. 

30. 

39 

0 

40 

1. 

0 

1. 

42 

-30. 

1200. 

30. 

39 

0 

40 

-1. 

0 

1. 

43 

-30. 

1200. 

-30. 

39 

0 

40 

-1. 

0 

-1. 

44 

30. 

1200. 

-30. 

39 

0 

40 

1. 

0 

-1. 

45 

30. 

1248. 

30. 

41 

30 

42 

1. 

1. 

1. 

46 

-30. 

1248. 

30. 

41 

30 

42 

-1. 

1. 

1. 

47 

-30. 

1248. 

-30. 

41 

30 

42 

-1. 

1. 

-1. 

48 

30. 

1248. 

-30. 

41 

30 

42 

1. 

1. 

-1. 

49 

90. 

1248. 

0. 

0 

31 

32 

1. 

1. 

0 

50 

-90. 

1248. 

0. 

0 

31 

32 

-1. 

1. 

c 

51 

60. 

1224. 

15. 

33 

34 

35 

1. 

1. 

1. 

52 

-60. 

1224. 

15. 

33 

34 

35 

-1. 

1. 

1. 

53 

60. 

1248. 

15. 

36 

B 

38 

1. 

1. 

1. 

54 

-60. 

1248. 

15. 

36 

If 

-1. 

1. 

1. 

55 

—60  . 

1224. 

-15. 

33 

34 

-1. 

X. 

-1. 

56 

60. 

1224. 

-15. 

33 

34 

35 

1. 

1. 

-1. 

12 

-60. 

60. 

1248. 

1248. 

;  -15. 

-15. 

It 

3? 

1! 

-1. 

1. 

1. 

1. 

-1. 

-1. 

30000000.  .1 

-15000 

.  20000. 

4. 

1 

1  5 

1  1 

10. 

2 

2  6 

1  1 

10. 

3 

3  7 

I  1 

10. 

4 

4  8 

1  1 

10. 

5 

5  9 

1  2 

10. 

6 

6  10 

1  2 

10. 

t 

7 

7  11 

1  2 

10. 

. 

8 

3  12 

1  2 

10. 

9 

9  13 

1  3 

10. 

.. _ 

H 

H  tt 

t  1 

M: 

67 


Table  XVIII  (cont'd) 


12 

12 

16 

1 

3 

10. 

13 

13 

17 

1 

4 

10. 

14 

14 

16 

1 

4 

10. 

15 

15 

19 

1 

4 

10. 

16 

16 

20 

1 

4 

10. 

17 

17 

21 

1 

5 

10. 

18 

18 

22 

1 

5 

10. 

19 

19 

23 

1 

5 

10. 

20 

20 

24 

1 

5 

10. 

21 

21 

25 

1 

6 

10. 

22 

22 

26 

1 

6 

10. 

23 

23 

27 

1 

6 

10. 

24 

24 

28 

1 

6 

10. 

25 

25 

29 

1 

7 

10. 

26 

26 

30 

1 

7 

10. 

27 

27 

31 

7 

10. 

28 

28 

32 

7 

10. 

29 

29 

33 

8 

10. 

30 

30 

34 

I 

8 

10. 

31 

31 

35 

8 

10. 

32 

32 

36 

8 

10. 

33 

33 

37 

9 

10. 

34 

34 

38 

9 

10. 

35 

35 

39 

9 

10. 

36 

36 

40 

9 

10. 

37 

37 

41 

X 

10 

10. 

38 

38 

42 

. 

10 

10. 

39 

39 

43 

10 

10. 

40 

40 

44 

'[ 

10 

io. 

41 

41 

45 

11 

10. 

42 

42 

46 

11 

10. 

43 

43 

47 

11 

10. 

44 

44 

48 

l 

Ll 

10. 

45 

5 

6 

l 

12 

10. 

46 

6 

7 

l 

12 

10. 

47 

7 

a 

l 

12 

10. 

48 

a 

5 

i 

12 

10. 

49 

9 

10 

l 

13 

10. 

50 

10 

ll 

i 

13 

10. 

51 

11 

12 

l 

13 

10. 

52 

9 

12 

l 

13 

10. 

53 

13 

14 

l 

14 

10. 

54 

14 

15 

l 

14 

10. 

55 

15 

16 

i 

14 

10. 

56 

13 

16 

l 

14 

10. 

57 

17 

18 

i 

15 

10. 

58 

18 

19 

l 

15 

10. 

59 

19 

20 

i 

15 

10. 

60 

17 

20 

l 

15 

10. 

61 

21 

22 

l 

16 

10. 

62 

22 

23 

l 

16 

10. 

63 

23 

24 

l 

16 

10. 

64 

21 

24 

l 

16 

10. 

65 

25 

26 

i 

17 

10. 

66 

26 

27 

l 

17 

10. 

67 

27 

28 

l 

17 

10. 

68 

25 

Is 

Z 

17 

10. 

69 

29 

30 

l 

18 

10. 

70 

30 

31 

l 

18 

10. 

71 

31 

32 

l 

18 

10. 

72 

29 

'  32 

l 

18 

:  0. 

73 

33 

34 

i 

19 

10. 

74 

34 

35 

i 

19 

10. 

75 

35 

36 

l 

19 

10. 

76 

33 

36 

i 

19 

10. 

77 

37 

38 

l 

20 

10. 

78 

38 

39 

l 

20 

io. 

79 

39 

40 

l 

20 

10. 

80 

37 

40 

l 

20 

10. 

81 

41 

42 

l 

21 

10. 

82 

42 

43 

l 

22 

10. 

83 

43 

44 

l 

21 

10. 

84 

41 

44 

i 

22 

10. 

85 

45 

46 

l 

23 

10. 
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Table  XVIII  (cont'd) 


86 

87 

83 

89 

90 

91 

92 

93 

94 

95 

96 

97 

98 

99 
100 
101 
102 

103 

104 

105 

106 

107 

108 
109 

M 

Hi 

tt* 

H  7 

119 

120 
121 
122 
123 

125 

& 

iff 

130 

131 

132 

133 

134 

135 

136 

137 

138 

139 

140 
14  1 

142 

143 

144 

145 

46 

47 

\t% 

li? 

ill 

'll 

.57 


46 

47 

1 

24 

10. 

47 

48 

1 

23 

10. 

45 

48 

1 

24 

10. 

1 

6 

1 

25 

10. 

2 

5 

1 

25 

10. 

2 

7 

1 

25 

10. 

3 

6 

1 

25 

1  0. 

3 

3 

I 

25 

10. 

4 

7 

1 

25 

10. 

4 

5 

1 

25 

10. 

1 

8 

1 

25 

10. 

5 

10 

1 

26 

10. 

6 

9 

1 

26 

10. 

6 

11 

1 

26 

10. 

7 

10 

1 

26 

10. 

7 

12 

1 

26 

10. 

8 

11 

1 

26 

10. 

8 

9 

1 

26 

10. 

5 

12 

l 

26 

10. 

9 

14 

1 

27 

10. 

10 

13 

1 

27 

10. 

10 

15 

1 

27 

10. 

11 

14 

1 

27 

10. 

it 

16 

1 

27 

1  0. 

12 

15 

l 

27 

10. 

12 

13 

1 

27 

10. 

9 

16 

1 

27 

10. 

13 

18 

1 

28 

10. 

14 

17 

i 

28 

10. 

14 

19 

i 

28 

10. 

15 

18 

i 

28 

10. 

15 

20 

i 

28 

10. 

16 

19 

i 

28 

10. 

16 

17 

i 

28 

10. 

13 

20 

i 

28 

10. 

17 

22 

i 

29 

10. 

18 

21 

i 

29 

10. 

18 

23 

i 

29 

10. 

19 

22 

i 

29 

10. 

19 

24 

i 

29 

io. 

20 

23 

i 

2? 

10. 

1% 


20 

17 

21 

22 

23 

ll 

24 
21 

\l 

26 

27 

27 

28 
28 

25 

fo 

1? 

31 

|| 

II 

if 

ll 

36 


21 

24 

25 

n 

ii 

28 

30 

29 

31 

30 

32 

31 
29 

32 

I 

34 

35 

u 

I! 

B 

i? 


10. 

10. 

10. 

10. 

10. 

10. 

10. 

10. 

10. 

10. 

10. 

10. 

10. 

10. 

10. 

10. 

10. 

10. 

10. 

0. 

0. 

0. 

0. 

10. 

10. 

10. 

10. 

10. 

io. 

10. 
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Table  XVIII  (cont'd) 


160 

161 

162 

«63 

164 

165 

166 

167 

168 

169 

170 

171 

172 

173 

174 

175 

176 

177 

178 

179 
ISO 
181 
132 

183 

184 

185 

186 

187 

188 

189 

190 

191 

192 

193 

194 

195 

196 

197 

198 

199 

200 
201 
202 

203 

204 

205 

206 

207 

208 

209 

210 
211 

hi 

f  15 
216 
217 


33 

37 

38 

38 

39 

39 

40 

40 
37 

41 

42 

42 

43 

43 

44 
44 
41 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 
50 
49 

49 

50 

50 
49 

45 

46 

47 

48 

51 

52 

55 

56 

53 

54 

51 

52 

53 
47 

45 

46 

43 

41 

42 

44 

5 

6 
9 


218 

14 

219 

17 

220 

13 

221 

21 

222 

2? 

223 

75 

224 

26 

225 

29 

226 

30 

22  r 

33 

228 

34 

229 

37 

230 

38 

41 

HI 

£ 

40 

1 

33 

42 

1 

34 

41 

1 

34 

43 

1 

34 

42 

1 

34 

44 

1 

34 

43 

1 

34 

41 

1 

34 

44 

1 

34 

46 

1 

35 

45 

1 

35 

47 

1 

35 

46 

1 

35 

48 

1 

35 

47 

1 

35 

45 

1 

35 

48 

1 

35 

51 

1 

36 

52 

1 

36 

55 

1 

36 

56 

1 

36 

53 

1 

37 

54 

1 

37 

57 

1 

37 

58 

1 

37 

51 

1 

38 

52 

1 

38 

55 

1 

38 

56 

1 

38 

53 

1 

39 

54 

1 

39 

57 

1 

39 

58 

1 

39 

51 

1 

40 

52 

1 

40 

55 

1 

40 

56 

1 

40 

53 

1 

41 

54 

1 

41 

57 

1 

41 

58 

l 

41 

58 

1 

42 

57 

1 

42 

56 

1 

43 

55 

L 

43 

48 

l 

44 

54 

l 

44 

58 

1 

44 

57 

l 

44 

52 

l 

45 

5o 

1 

45 

55  i 

L 

45 

51 

1 

45 

7 

1 

46 

8  ] 

L 

46 

11 

£ 

47 

12 

i 

47 

15 

i 

48 

16 

i 

48 

19  ; 

f 

49 

20  ] 

l 

49 

23  1 

l 

50 

24  1 

L 

50 

27  i 

£ 

51 

28  1 

i 

51 

31  ; 

£ 

52 

32  ] 

i 

52 

34 

t 

53 

36 

53 

39 

£ 

54 

40 

54 

43 

£ 

55 

44 

i 

55 

47 

i 

56 

10. 

10. 

10. 

10. 

10. 

10. 

10. 

10. 

10. 

10. 

10. 

10. 

10. 

10. 

10. 

10. 

10. 

10. 

10. 

10. 

10. 

10. 

10. 

10. 

10. 

10. 

10. 

10. 

10. 

10. 

10. 

10. 

10. 

10. 

10. 

10. 

10. 

10. 

10. 

10. 

10. 

10. 

10. 

10. 

10. 

■  0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
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'able  XVIII  (cont'd) 


234 
1 
2 

3 

4 
2 

49 

50 
2 

49  6000. 

50  -6000. 
2 


46 

l 

1 

1 

1 

6000. 

6000. 


48 

1 

1 

1 

1 


1 
1 
1 
1 
1 

-20000. 

-20000. 

-20000. 

-20000. 


56 


10. 


0.0 

0.0 

0.0 

0.0 


5888:  -!8888:  hi 


49 

50 
2 

49  3000. 

50  3000. 

2 

49  -3000. 

50  -3000. 

49  200. 

50  200. 

25.  25. 

25.  25. 

25.  25. 

25.  25. 

25.  25. 

25.  25. 

25.  25. 

0.000001  .000001 
0.000001  .000001 
0.000001  .000001 
0.000001  .000001 
.000001  .000001 
.000001  .000001 


-10000. 

-10000. 

10000. 

10000. 


10. 


-5000. 

5000. 

5  000. 
-5000. 


25. 

25. 

25. 

25. 

25. 

25. 

25. 

.000001 

.000001 

.000001 

.000001 

.000001 

.000001 

.000001 

100. 

100. 

100. 

100. 


108: 
100. 
111. 
93. 
720. 
39. 
60. 

0. 

0. 

0. 

0. 

0. 

2000. 

2000. 

2000. 

2000. 

2000. 


if: 

II: 

25. 

25. 

25. 

.003001 
.000001 
000001 
000001 


. 000001 
.000001 
100. 
100. 
100. 
100. 
ICO. 
100. 
100. 
120. 

84. 

66. 

1080. 

15. 

0. 

0. 

0. 

0. 

0. 

2000. 

2000. 

2000. 

1888: 


25. 

25. 

25. 

25. 

25. 

25. 

25. 

.000001 

.000001 

.000001 

.OoW 

.000001 

.000001 

100. 

100. 

100. 

100. 

M: 

100. 

US: 

11: 

1248. 

3. 

0. 

0. 

0. 

0. 

2030. 

2000. 

2000. 

1888: 


25. 

25. 

25. 

25. 

25. 

25. 

25. 

.000001 

.000001 

.000001 

.000001 

.000001 

.000001 

.000001 

100. 

100. 

100. 

100. 

100. 

100. 

100. 

102. 

84. 

340. 

1248. 

60. 

0. 

0. 

0. 

0. 

0. 

2000. 
2000. 

?00C 


25. 

25. 

25. 

25. 

25. 

II: 

.000001 

.000001 


.  00000, 
000001 
.000001 
.000001 
100. 

100. 

100. 

100. 

100. 

100. 

100. 

»!• 

fl;«. 

30. 

0. 

0. 

0. 

0. 

0. 


1888 

200“ 

188 


25. 

25. 

25. 

25. 

25. 

25. 

25. 

.000001 
-OOOCCi 
.000001 
.000001 
.000001 
.000001 
. 000001 
100. 
ICO. 
100. 
100. 
ICO. 
100. 
100. 
102. 
600. 

46. 

90. 

30. 
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APPENDIX  3 


PROGRAM  ORGANIZATION 


A.  DESCRIPTION 

The  program  organization  is  laid  out  in  the  following 
flow  charts .  The  main  driver  program  is  arranged  to  call  a 
subdriver,  XMSADT,  and  the  optimizer  of  the  user's  choice. 
All  changes  for  replacing  an  optimizer  occur  in  MSADT.  This 
allows  for  easy  testing  of  several  optimizers  on  the  same 
problem. 

XMSADT  may  be  called  from  the  main  for  input,  analysis, 
and  output.  Printed  output  may  vary  as  the  user  requires. 

A  complete  listing  of  all  subroutines  and  their  functions  is 
given  in  Table 
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Table  XIX 


ANALYSIS  CONTROL 

ANALYSIS  ICALC=2 


ICALC=2 


Table  XX 


OUTPUT  CONTROL 
OUTPUT  FOR  INFO=0/l/2 

,  I 

Ij 


ICALC=  3 


(PRINT) 
-GRADIENTS 
SPECIFIED  BY 
INPUT  CONTROL 
(IDVCLC) 


( 


Table  XXI 

PROGRAM  FLOW  CHART 


..  _LIM£JJT)_  _  _  _ 

-  SADT46 


SADTO 1- 


SADT02 


SADT11 


I 


SADT04- 
SADT07 
SADTO 8 


SADTO 3 
SADTO 5 


SADT1-5—  SADT16 
\ — SADT23 


“SADT06  —  SADT19 
1  SADT6 1”  "SADT4  8  ' 

-SADTO 5 

-  SADT5  0  _ SADTO  8 

I-  SADT4  9  - [~ 

I - EIGRQ1 

I - - 1 


EIGR02  EIGR03 


(ANALYSIS ) 


L 

r 


‘SADT7  3“ 

SADT75 

SADT76 

SADT77 


SADT36 

SADT20' 

SADT21 


(GRADIENTS) 


C-" 

C= 


SADT74 

mm 

SADT79- 

SADT17 

SADT62 

SADT71 


SADTO  3 
SADT4  8 

SADT37 


(OUTPUT) 


TABLE  X XII 


FILENAME 

EIGBO 1 
EIGR02 
EIGB03 
SADTO  1 

SADT02 

SADT03 

SADTO  4 

SADT05 

SADT06 

SADT07 

SADT08 

SADT09 

SADT10 

SADT11 

SADT14 

S A DTI  5 

S A DTI 6 
SADT17 
SADT18 
SADT19 


SOBBOOTINE  OSES 


TITLE 


SOLVES  EIGEN  7AL0E  PROBLEM  |A-ALAMBD*B|  X=0 

SOLVES  EIGENVALUE  PROBLEM 

SOLVES  EIGEN VALOE  PROBLEM 

THIS  ROOTINE  READS  AND  PRINTS  INPOT  DATA 
AND  ORGANIZE S  PS EUDO-DY NAMIC  STORAGE 
ALLOCATION 

BOILDS  VECTORS  JC  AND  UK  FOR  FINITE  ELEMENT 
STROCTORAL  ANALXSIS 

BOILDS  THE  ELEMENT  STIFFNESS  MATRIX  FOR 
TROSS  FINITE  ELEMENT  NO.  II 

BOILDS  ELEMENT  STIFFNESS  MATRICES  AND  ADDS 
THEM  TO  THE  GLOBAL  STIPFNESS  MATRIX 

SUPERIMPOSES  THE  ELEMENT  STIFFNESS  VECTOR  EK 
OF  ELEMENT  1 1  ON  THE  GLOBAL  VECTOR  AK 

BUILDS  GLOBAL  LUMPED  MASS  MATRIX 

LO  DECOMPOSES  SYMMETRIC. POSITIVE-DEFINITE. 
SPARCE  MATRICES, THE  UPPER  TRIANGLE  OF  WHitH 
IS  STORED  IN  VECTOR  AK  WITH  LEADING  ZEROES 
OMITTED 

FORWARD  AND  BACK  SUBSTITUTES  TO  YIELD  A  SOLN 
TO  A  SET  OF  LINEAR  EQNS ( DECOMPOSED  W/ 

SADT07) 

PRINTS  ALL  JOINT  DISPLACEMENTS  FOR  EACH 
LOAD  CONDITION  OF  A  FINITE  ELEM.  STRUCTURE 

PERFORMS  TRUSS  FIXED-GEOMETRY  DESIGN 

ROUTINE  TO  ORGANIZE  ANALYSIS 

CALCULATES  PARTIAL  DERIVATIVE  OF  TRUSS  ELEM 
K  MATRIX  WRT  COORD  DOF  II  AT  NODE  J1 

CALCULATESS  VALUES  OF  ALL  DESIGN  AND  BEHAV- 
IORIAL  CONSTRAINTS  AS  DEFINED  BY  PROGRAM 
"SAD" 

CALCULATES  STRESS  IN  BAR  ELEMENT  II  UNDER 
LOAD  CONDITION  JJ 

PRINTS  FORCES  OR  STRESSES  IN  BAR  ELEMENTS 


CALCULATES  GRADIENT  INFO  IN  FINITE  ELEM 
STRUCTURAL  ANALYSIS 

ADDS  ELEHEMENT  MASS  AA  OF  ELEMEMENT  II 
TO  GLOBAL  NASS  MATRIX 
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SADT20 

SADT21 

SADT22 

SADT23 

SADT26 

SADT31 

SADT34 

SADT35 

SADT36 

SADT37 

SADT40 

SADT41 

SADT42 

SADT43 

SADT44 

SADT45 

SA0T46 

SADT48 

SADT49 

SADT50 

SADT51 


CALCULATES  THE  PARTIAL  DERIVATIVES  OF  THE 
DISPLACEMENTS  WRT  INDEPENDENT  DESIGN 
VARIABLE  II 


CALCULATES  STRESS  AND  GRADIENT  OF  STRESS 
IN  BAR  ELEM  IllJNDER  LOAD  COND  JJ  GIVEN 
DISPLACEMENT  VECTOR  PL  AND  GRADIENTS  OF 
THE  DISPLACEMENTS  IN  VECTOR  DU 


PERFORMS  STRESS-RATIO  DESIGN  OF  A  FIXED 
GEOMETRY  TRUSS 

CRATES  WEIGHT  OF  TRUSS  GIVEN  AREAS  AND 

CALCULATES  THE  GRADIENT  OF  WEIGHT  WRT 
RECIPROCAL  DESIGN  VARIABLES 

CALCULATES  GRADIENT  INFORMATION  AND 
STORES  IN  ARRAY  TAY 

CONMIN  EXTERNAL  FOR  TRUSS  FIXED  GEOMETRY 
DESIGN 

CALCULATES  GRADIENT  INFORMATION  AND 
STORES  IN  ARRAY  TAY 

CALCULATES  t  XEIG-T*AM* XEIGI  FOR 
UCULATIONS  IN  FR 


GRADIENT  CALCULATIONS 
CONSTRAINTS 


EQUENCY 


CALCULATES  EIGENVALUE  GRADIENT  INFORMATION 
IN  PINITE  ELEMENT  STRUCTURAL  ANALYSIS  AND 
DESIGN 

CALCULATES  GRADIENTS  OF  DISPLACEMENTS 
WRT  COORDINATE  DESIGN  VARIABLE  II 

CALCULATES  GRADIENT  INFORMATION  WRT 
COORDINATE  VARIABLES 

CALCULATES  GRADIENT  OF  WEIGHT  WRT 
COORDINATE  DESIGN  VARIABLES 

CALCULATES  GRADIENTS  OF  OBJECTIVE  AND 
ACTIVE  CONSTRAINTS  IN  OPTIMUM  GEOMETRY 
DIRECTION  FINDING  PROBLEM 

PERFORMS  TRUSS  FIXED  GEOMETRY  DESIGN  AND 
EVALUATES  CONSTRAINTS 

PERFORMS  TRUSS  OPTIMUM  GEOMETRY  DESIGN 

READS  INPUT  INFORMATION  FOR  BAR  ELEMENTS 

BUILDS  THE  ELEMENT  MASS  MATRIX  FOR  A 
TRUSS  FINITE  ELEMENT 

SOLVES  REAL  EIGENVALUE  PROBLEMS  USING 
THE  SUBSPACE  ITERATION  METHOD 


BUILDS  INITIAL  SET  OF  BASIS  VECTORS 
FOR  El  GEN  SOL  UTION  BY  REDUCED  BASIS 
METHOD 

CALCULATES  STRESS  AND  GRADIENT  OF  STRESS 
IN  BAR  ELEMENT  II  UNDER  LOAD  CONDITION 
JJ,  GIVEN  DI SLPL ACEH ENTS  CONTAINED  IN 
VECTOR  PL  AND  GRADIENT  OF  D7C PLACEMENTS 
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i 


IN  VECTOR  DO 


SADT52 

SADT53 

SADT54 

SADT6 1 
SADT62 

SADT71 

SADT73 

SADT74 

SADT75 

SADT76 

SADT77 

SADT79 

MSADT 

XHSADT 


S?§I8?B!iJ-DoiI!IiI8?ibSE*aca  IS 

PRINTS  MEMBER  INFORMATION  FOR  BAR  ELEM 

CALCULATES  AND  PRINTS  CENTER-OF-GR A VITY 
AND  INERTIAL  PROPERTIES  OF  FINITE  ELEM 
STRUCTURE 

BUILDS  GLOBAL  MASS  MATRIX 

PRINTS  NEIG  EIGENVALUES  STORED  IN 
EIGV AL ,  AND  THEIR  CORRESPONDING 
EIGENVECTORS  STORED  IN  XEIG 

PRINTS  G  VECTOR  OF  CONSTRAINTS 

CALCULATES  GRADIENT  INFORMATION 

PRINTS  TAY  ARRAY  OF  GRADIENTS 

CALCULATES  GRADIENT  OF  HEIGHT  WITH 
RESPECT  TO  AREA  DESIGN  VARIABLES 

CALCULATES  GRADIENT  OF  HEIGHT  HITH 
RESPECT  TO  COORDINATE  DESIGN  VARIABLES 

CALCULATES  GRADIENT  INFORMATION 

PRINTS  COORDINATE  INFORMATION 

DRIVER  PROGRAM  FOR  USING  THE  ABOVE 
SUBROUTINES.  MSADT  MAY  BE  COUPLED 
TO  OPTIMIZER  OF  USER'S  CHOICE. 

SUBDRIVER  PROGRAM  FOR  COUPLING 
SADT  ROUTINES  TO  MSADT 


LIST  OF  REFERENCES 


1.  Vanderplaats ,  G.  N. ,  "CONMIN  -  A  Fortran  Program  for 
Constrained  Function  Minimization  -  User's  Manual," 

NASA  TMX-62,  282,  August  1973. 

2.  Felix,  J.  and  Vanderplaats,  G.  N.,  Configuration 
Optimization  of  Trusses  Subject  to  Strength,  Displace¬ 
ment  and  Frequency  Constraints,  proc.  ASME  Computer 
Engineering  Conference  and  Exhibit ,  San  Diego,  August 
15-19,  1982. 

3.  Arora,  J.  S.  and  Haug,  E.J.,  "Methods  of  Design  Sensi¬ 
tivity  Analysis  in  Structural  Optimization , 11  AIAA  Journal, 
Vol.  17,  No.  9,  September  1979,  pp .  970-174. 

4.  Vanderplaats,  G.  N.,  "Comment  on  Method  of  Design 
Sensitivity  Analysis  in  Structural  Optimization,"  AIAA 
Journal,  Vol.  18,  No.  11,  November  1976  ,  pp .  1406-1407 . 

5.  Nelson,  R.  B.,  "Simplified  Calculation  of  Eigenvector 
Derivatives,"  AIAA  Journal,  Vol.  14,  No.  9,  September 
1976,  pp.  1201-1505. 


80 


INITIAL  DISTRIBUTION  LIST 


No 


1.  Defense  Technical  Information  Center 
Cameron  Station 

Alexandria,  Virginia  22314 

2.  Library,  Code  0142 
Naval  Postgraduate  School 
Monterey,  California  93940 

3.  Professor  G.  N.  Vanderplaats ,  Code  69Vn 
Department  of  Mechanical  Engineering 
Naval  Postgraduate  School 

Monterey,  California  93940 

4.  Professor  Young  Shin,  Code  69 
Department  of  Mechanical  Engineering 
Naval  Postgraduate  School 
Monterey,  California  93940 

5.  Dr.  Jim  Bennett 

GM  Research  Laboratory 
Warren,  Michigan  48090 

6.  Dr.  Richard  L.  Bennett 
Bell  Helicopter  Textron 
P.0.  Box  482 

Forth  Worth,  Texas  76101 

7.  Dr.  Richard  H.  Gallegner 

Dean  of  College  of  Engineering 
University  of  Arizona 
Tucson,  Arizona  85721 

8.  Dr.  Erwin  Johnson 

Northrop  Corporation  (Aircraft  Group) 

1  Northrop  Avenue 
Hawthorne,  California  90250 

9.  Dr.  Hirokazu  Miura 

NASA  Ames  Research  Center 
MS  237-11 

Moffett  Field,  California  94035 


Copies 

2 

2 

10 

1 

1 

1 

1 

1 

1 


31 


No.  Copies 
1 


10.  Dr.  Lucien  A.  Schmit 
Department  of  Mechanical  5  Structural  Engineering 
6731  3oelter  Hall 

University  of  California  Los  Angeles 
Los  Angeles,  California  90024 

11.  Dr.  Jarek  Sobieski  1 

NASA  Langley  Research  Center 

MS  24  3 

Hampton,  Virginia  23665 

12.  LT  James  E.  "itzgerald  2 

122  Pennsylvania  Avenue 

Lebanon,  Tennessee  37Q82 

13.  Chairman,  Code  69  1 

Department  of  Mechanical  Engineering 

Naval  Postgraduate  School 
Monterey,  California  93940 


1  i 

i  ' 
I  1 


I 


j 

i 


i 

i 

j 


32 


